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THE  INFLUENCE  OF  CONVECTION  VELOCITY  ON  THE 
TURBULENT  WALL  PRESSURE  WAVENUMBER-FREQUENCY  SPECTRUM 


INTRODUCTION 

Empirically  derived  models  of  the  turbulent  wall  pressure  wavenumber-frequency  spectrum  are 
required  in  problems  of  flow-induced  vibration  and  noise.  Of  particular  interest  is  the  energy  in  the 
subconvective  and  low-wavenumber  region,  which  excites  the  lower  order  vibrational  modes  of 
structures.  A  direct  measurement  of  C>(ki,k2,(jL))  requires  an  array  of  many  closely  spaced  sensors 
from  which  data  are  acquired  simultaneously  in  time.  Due  to  the  extensive  number  of  sensors  and 
signal  processing  required,  very  few  such  measurements  have  been  made  in  experimental  flow 
facilities.  The  inherent  problem  of  spatial  aliasing  due  to  sensor  spacing  limits  the  wavenumber 
and  frequency  bandwidths  of  the  measurements,  as  discussed  by  Keith  and  Abraham  (1994). 
Models  of  O(ki,k2,co)  are  therefore  often  derived  from  measurements  of  the  cross  spectrum  made 
with  a  relatively  small  number  of  sensors.  Uncertainties  result  from  the  limited  number  of  spatial 
separations  as  well  as  from  the  size  of  the  pressure  sensors,  which  determines  the  separations  that 
are  possible. 

A  convection  or  phase  velocity,  uc,  may  be  determined  from  the  streamwise  cross-spectrum 
phase  (t)(^,co)  by  use  of  the  relationship  (1)(^,«))  =  co^Uc-  It  is  commonly  assumed  in  estimating  the 
wavenumber-frequency  spectrum  that  uc  is  independent  of  In  fact,  Uc  varies  with  ^  due  to  the 
variation  in  convection  velocities  and  spatial  decay  rates  of  the  pressure-producing  turbulent 
structures  in  the  inner  and  outer  regions  of  a  turbulent  boundary  layer.  Here,  we  first  formulate  the 
problem  in  the  spatial  domain  using  the  Corcos  (1963)  cross-spectrum  model.  We  then  compare 
experimental  and  numerical  data  for  uc  from  investigations  covering  a  wide  range  of  Reynolds 
numbers  in  order  to  determine  effective  scaling  parameters.  The  wavenumber-frequency  spectra 
are  next  estimated  by  transforming  the  cross  spectrum,  taking  into  account  the  dependence  of  uc 
upon  Comparisons  with  recent  numerical  results  are  shown  to  be  favorable,  and  the  effects  on 
the  attenuations  in  the  autospectra  due  to  finite  sensor  size  are  determined. 

MODELING  THE  WALL  PRESSURE  CROSS  SPECTRUM 

The  streamwise  wall  pressure  cross  spectrum  may  be  expressed  as 

G(^,cd)  =  IG(^,(o)l  e-i<t>(^>“) .  (1) 
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For  a  fixed  frequency  coq,  assuming  uc  is  independent  of  ^  and  IG(^,C0o)l  decays  exponentially 
with  equation  (1)  becomes 


G(^,coo)  =  ’ 

where  a  is  a  decay  constant.  The  spectrum  O(ki,C0o)  is  then 

O(ki,coo)  =  <J>(coo)/27t  r  e(«l^l-i(kl  +  ®o/uc(«o))Od^  . 


(2) 


(3) 


Using  MACSYMA  (Rand  1984),  we  obtain  the  closed-form  solution 

<D(ki,coo)  =  (d>(03o)/2tt)(-2a/(a2  -i-  (ki  kc)^)) .  (4) 

If  Uc  is  allowed  to  vary  with  ^  as  well  as  coq  in  equation  (3),  a  spatial  modulation  of  the  phase  term 
results,  and  in  general  the  integral  cannot  be  evaluated  in  closed  form.  Assuming  uc  to  be 
independent  of  we  define  a  =  y  CDo/uc(<»o).  and  express  equation  (2)  in  terms  of  the  similarity 
variable  =  O)o^uc(coo)  such  that 

G(^*,coo)  =  4)(coo)e(^^*'‘^  ■  (5) 

The  Corcos  (1963)  cross-spectrum  model  is  separable  in  the  streamwise  and  spanwise  coordinates. 
Equation  (5)  represents  the  streamwise  component  of  the  model.  Recently  Farabee  and  Casarella 
(1991)  and  Keith  and  Barclay  (1993)  have  shown  that  the  similarity  scaling  holds  over  a  limited 
frequency  range.  Although  Corcos  did  not  transform  his  cross-spectrum  model  to  the 
wavenumber  domain,  equation  (4)  with  a  replaced  by  ykc  is  commonly  referred  to  as  the  Corcos 
model  of  <I)(ki,C0o)-  Corcos  defined  a  convection  velocity  in  terms  of  the  phase  of  the  cross 
spectmm,  and  stated  that  "uc  is  found  to  be  a  function  of  co  and  to  be  almost  independent  of  ^  and 
T]."  This  conclusion  was  likely  based  upon  the  limited  amount  of  convection  velocity  data  available 
at  that  time. 

Farabee  and  Casarella  (1991)  showed  that  the  similarity  scaling  proposed  by  Corcos  was 
effective  for  collapsing  their  cross-spectral  data,  taking  into  account  the  variation  of  Uc  with 
However,  the  influence  of  the  variation  of  uc  with  ^  on  the  wavenumber-frequency  spectrum  has 
not  been  investigated  to  date.  A  criticism  of  the  model  for  <I>(ki,coo)  given  by  equation  (4)  is  an 
apparent  overprediction  of  spectral  levels  at  subconvective  and  lower  wavenumbers.  Here,  higher 
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spectral  levels  in  these  regions  will  be  shown  to  result  from  constraining  uc  to  be  constant  when 
transforming  the  cross  spectrum.  Taking  into  account  the  variation  of  uc  with  ^  will  be  shown  to 
lead  to  improved  estimates  of  <I>(ki,0)o). 


COMPARISON  OF  WALL  PRESSURE  CONVECTION  VELOCITIES 

The  dependence  of  the  convection  velocity  on  the  spatial  separation  and  effective  scaling  laws 
will  be  determined  by  comparisons  of  previous  measurements.  The  dominant  trend  in  the  data  is 
an  increase  in  uc(^  ,co)  as  ^  increases  and  to  decreases.  This  results  from  the  locations  of  the  waU- 
pressure-producing  turbulent  structures  in  the  boundary  layer  as  well  as  from  their  rate  of  spatial 
decay.  The  larger  stmctures  at  locations  farther  from  the  wall  convect  faster  and  are  coherent  over 
larger  distances  than  those  near  the  wall.  For  larger  separations,  the  contributions  from  the  larger 
structures  therefore  dominate  uc(^,oo).  At  smaller  separations,  the  contributions  from  the  smaller 
structures  coherent  over  those  distances  cause  a  reduction  in  uc(^,Cl)).  The  increase  in  uc(^,co)  with 
decreasing  co  is  also  consistent  with  these  effects. 

Here  we  focus  on  the  spatial  variation  of  uc  for  discrete,  narrowband,  and  broadband  frequency 
ranges.  We  compare  data  from  six  investigations  that  cover  two  orders  of  magnitude  in  R0.  The 
data  of  Farabee  and  Casarella  (1991)  were  determined  from  the  phase  of  the  cross  spectrum  at 
discrete  frequencies;  those  of  Willmarth  and  Wooldridge  (1962)  and  Keith  et  al.  (1991)  were 
determined  from  narrowband  correlation  measurements;  and  those  of  Choi  and  Moin  (1990), 
Schloemer  (1967),  and  Bull  (1967)  from  broadband  correlation  functions.  Farabee  and  Casarella 
(1991)  determined  scaling  laws  for  wall  pressure  autospectra  and  the  regions  of  the  boundary  layer 
that  contribute  at  particular  frequency  ranges.  Use  of  these  length  and  time  scales  results  in 
scalings  given  by  uc/Uq  vs.  ^5*,  uc/ux  vs.  ^5,  and  uc/ux  vs.  ^ux/v. 

For  comparing  the  data  using  the  outer  scaling  uc/Uo  vs.  ^5*,  we  consider  the  frequency  bands 
scaled  as  cl)5*AJo-  The  higher  discrete  frequency  of  Farabee  and  Casarella  falls  at  the  center  of  the 
narrow  bandwidth  of  Keith  et  al.  (at  both  R0).  These  data  display  a  reasonable  collapse,  as  shown 
in  figure  1.  Although  the  narrow  bandwidth  of  Willmarth  and  Wooldridge  lies  above  the  higher 
discrete  frequency  of  Farabee  and  Casarella,  these  data  coUapse  well.  The  numerical  simulation  of 
Choi  and  Moin  covers  approximately  the  same  broad  band  as  Bull's.  Choi  and  Moin's  data  are 
above  Bull's  at  the  lower  separations,  but  converge  at  the  higher.  The  broadband  data  of 
Schloemer  and  Bull  agree  well  where  Schloemer's  data  covers  a  slightly  higher  frequency  range 
than  do  Bull's.  The  low  discrete  frequency  of  Farabee  and  Casarella  lies  within  the  broad  band  of 
Bull.  The  two  data  sets  agree  well  for  the  larger  separations.  This  agreement  between  broadband 
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and  discrete  low-frequency  data  reflects  the  inherent  filtering  of  the  higher  frequency  contributions 
in  the  broadband  data  resulting  from  the  spatial  separations  of  the  sensors. 

The  outer  scaling  effectively  collapses  the  discrete,  narrowband  and  broadband  data  sets.  The 
dominant  trend  of  increasing  convection  velocity  with  increasing  separation  is  apparent.  We 
investigated  the  use  of  5  rather  than  6*  as  an  outer  length  scale,  and  found  no  significant  effect.  At 
low  values  of  ^/5*,  the  data  sets  at  higher  R0  tend  to  display  a  sharper  roll-off  than  Farabee  and 
Casarella’s  low-frequency  data.  For  the  scaling  uc/ux  vs.  ^/5  of  figure  2,  we  consider  the 
frequency  bands  scaled  as  (oS/ux-  Although  the  higher  discrete  frequency  data  of  Farabee  and 
Casarella  fall  at  the  center  of  the  narrow  bands  of  Keith  et  al.  at  both  Rq,  these  data  do  not  collapse 
well.  In  addition,  the  collapse  of  the  data  of  Keith  et  al.,  which  covers  the  same  narrow  band  at 
both  R0,  is  poor,  with  the  lower  R0  data  lying  below  the  higher.  Bull  and  Schloemer's  broad 
bands  are  comparable,  and  the  two  data  sets  collapse  fairly  well,  with  Schloemer's  data  slightly 
below  Bull's.  Choi  and  Moin's  broad  band  covers  the  lower  half  of  Bull's,  but  their  data  lies 
significantly  below  Bull's.  The  use  of  ux  as  a  velocity  scale  combined  with  an  outer  length  scale 
therefore  does  not  effectively  coUapse  data  over  a  wide  range  of  R0.  Data  at  lower  R0  are  mapped 
below  the  higher  R0  data,  resulting  in  the  trend  of  increasing  levels  with  increasing  R0. 

The  use  of  ux  as  a  velocity  scale  combined  with  the  inner  length  scale,  v/ux,  leads  to  an 
improved  collapse  of  the  data  of  figure  2,  as  shown  in  figure  3.  For  this  inner  scaling,  we 
consider  the  frequency  bands  scaled  as  cov/ux^.  The  low  discrete  frequency  data  of  Farabee  and 
Casarella  fall  at  the  center  of  the  narrow  band  of  Keith  et  al.  at  the  higher  R0.  These  data  collapse 
well.  The  high  discrete  frequency  data  of  Farabee  and  Casarella  falls  at  the  center  of  the  narrow 
band  of  Willmarth  and  Wooldridge.  The  Willmarth  and  Wooldridge  data  lie  above  those  of 
Farabee  and  Casarella.  Choi  and  Moin's  broadband  data  extend  to  significantly  higher 
frequencies,  cov/ux^,  than  the  other  data  sets.  Their  data  lies  significantly  below  those  of  the  other 
three  investigations.  Although  Schloemer's  data  covers  a  somewhat  higher  frequency  range  than 
do  Bull's,  their  data  collapse  fairly  well.  The  data  of  Farabee  and  Casarella  (low  co),  BuU,  Keith  et 
al.  (high  R0),  and  Willmarth  and  Wooldridge  collapse  well  at  larger  values  of  ^  uxA'. 

The  outer  variable  scaling  of  figure  1  collapses  all  of  the  data  sets  over  two  orders  of  magnitude 
in  R0  somewhat  more  effectively  than  does  the  inner  variable  scaling.  Based  upon  the  analysis  of 
Farabee  and  Casarella  (1991),  in  view  of  the  frequency  bands  as  given  in  table  1,  contributions  to 
the  wall  pressure  autospectra  for  these  cases  are  from  a  significant  or  entire  portion  of  the  log-law 
region.  The  effectiveness  of  outer  and  inner  variables  for  collapsing  the  convection  velocity  data  is 
consistent  with  the  sources  existing  in  the  log-law  region. 
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Figure  1.  Comparison  of  Convection  Velocities  Scaled  With  Outer  Variables 


Figure  2,  Comparison  of  Convection  Velocities  Scaled  With  Mixed  Variables 


Table  1.  Parameters  for  Comparison  of  Convection  Velocity  Data 


Investigation 

Re 

co8*/Uoll 

co6*/Uol2 

co6/ux!  1 

co6/uxl2 

cov/ux^ll 

cov/ux^l2 

Region 

Choi  &  Moin 
(1990),  R(^,t), 

285 

(DNS) 

180 

0.08 

3.55 

10.00 

450.00 

0.055 

2.484 

mid-high 

Farabee  & 

Casarella  (1991), 
O(^,C0) 

3,400 

(air) 

1,165 

0.24 

0.96 

40.8 

163.0 

— - 

0.035 

0.140 

— 

mid 

log-law 

Schloemer 
(1967).R(^.x) 
Broad  Band 

5,800 

(air) 

2,043 

0.115 

7.655 

20.772 

1,384.83 

0.010 

0.678 

mid-high 

BuU  (1967), 

Rlc  t) 

10,000 

(air) 

3,614 

0>038 

4.08 

8.24 

884.30 

0.002 

0.245 

mid-  log 
law 

Keith  et  al. 

(1991),  R(^,x) 

11.700 

21.700 
(water) 

4,182  1 

6,778 

0.47 

0.46 

1.41 

1.42 

100.09 

106.29 

300.28 

328.11 

0.024 
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Figure  3.  Comparison  of  Convection  Velocities  Scaled  With  Inner  Variables 
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ESTIMATED  WAVENUMBER-FREQUENCY  SPECTRUM  0*(Ki,Q) 


Defining  nondimensional  outer  variables  X  =  ^/5*,  Ki  =  -ki5*,  Uc  =  uc(C0o,^)/Uo,  O.  = 
coo5*AJo.  and  Kc  =  Q/Uc,  a  nondimensional  normalized  form  of  the  Corcos  model  is  given  by 

<D*(Ki,Q)  =  (-2y/Kc)(1/(72  +  (1  -  Ki/Kc)^)) .  (6) 

If  the  dependence  of  the  convection  velocity  on  spatial  separation  is  taken  into  account,  the 
spectrum  may  be  evaluated  as 

0*(Ki,Q)  =  J  e(Tli^X/Uc(X)|-i(-Ki  +  Q/Uc(X))X)dx  .  (7) 


The  sign  of  Ki  was  chosen  as  negative  (in  defining  nondimensional  variables)  to  place  the 
convective  ridge  at  positive  Ki.  A  value  for  the  decay  constant  y  of  -0.145  provided  an  accurate 
exponential  curve  fit  to  the  magnitude  of  Farabee  and  Casarella's  (1991)  cross  spectrum  for  Q  = 
0.24.  A  function  of  the  form  Uc(X)  =  aX^*  was  fit  to  the  low-frequency  data  of  Farabee  and 
Casarella  (1991),  resulting  in  a  standard  deviation  of  0.04.  For  X  >  81,  Uc(X)  was  taken  as 
constant  at  0.85.  This  curve-fit  results  in  a  sharper  roll-off  near  the  origin  than  that  supported  by 
Farabee  and  Casarella's  data,  and  more  accurately  models  the  broadband  data  of  Bull  and  the 
higher  frequency  data  of  Willmarth  and  Wooldridge,  as  shown  in  figure  4.  The  function  Uc(X)  = 
a(X-i-l)b  resulted  in  the  same  standard  deviation  of  0.04  for  Farabee  and  Casarella’s  low-frequency 
data,  and  gives  a  value  for  Uc(0)  of  0.645,  which  agrees  with  that  estimated  by  Farabee  and 
Casarella. 

Equation  (7)  was  evaluated  numerically,  with  Uc(X)  =  aX^*.  The  resulting  spectrum  of  figure  5 
exhibits  significantly  lower  levels  at  subconvective  and  lower  wavenumbers  and  increased  levels  at 
’  higher  wavenumbers,  in  comparison  to  the  Corcos  model  given  by  equation  (6)  with  a  constant 

value  of  0.85  for  Uc-  The  function  Uc(X)  =  a(X-i-l)t>  leads  to  the  same  trends  in  the  spectrum  as 
the  former,  with  the  effects  slightly  diminished.  For  comparison,  we  also  considered  the  function 
a(X+10)t'  with  less  of  a  roll-off  near  the  origin.  This  function  produced  the  same  trends  in  the 
spectrum,  with  the  effects  further  diminished.  A  fifth-order  polynomial  (not  shown  here)  was  also 
fit  to  the  data  of  Farabee  and  Casarella,  and  resulted  in  a  spectrum  that  agreed  quite  well  with  that 
obtained  using  the  function  a(X+l)b.  The  function  a(X-i-l)t’  was  found  to  most  accurately  model 
the  convection  velocity  data  and  provides  the  best  estimate  of  the  wavenumber-frequency  spectra 
for  this  case. 
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Figure  5.  Effects  of  the  Roll-Off  of  Uc(X)  for  Small  X  on  the 
Wavenumber-Frequency  Spectrum 


Comparison  of  the  cases  aX^  and  a(X+l)t>  provides  insight  into  the  sensitivity  of  the  spectrum 
to  small-scale  turbulence.  The  roll-off  of  Uc(X)  near  the  origin  results  from  small-scale  turbulent 
structures  with  lower  convection  velocities  and  rapid  spatial  decay.  For  X  >  1,  the  values  for 
Uc(X)  given  by  aX^  and  a(X-(-l)t’  agree  to  within  5  percent.  The  differences  in  the  roll-off  of 
Uc(X)  for  0  <  X  <  1  therefore  lead  to  the  differences  in  the  spectra  at  low  wavenumbers,  with  the 
sharper  roll-off  leading  to  lower  levels.  For  the  boundary  layer  of  Farabee  and  Casarella,  X  =  1 
corresponds  to  ^  =188v/ut.  A  decay  constant  of  0.145  implies  that  a  turbulent  pressure-producing 
structure  loses  84  percent  of  its  coherence  as  it  convects  and  decays  over  2X.  A  cut-off  may  then 
be  defined  as  X  <  lOOv/ux  for  turbulence  that  influences  Uc(X)  only  over  this  particular  spatial 
range. 

At  the  low  frequency  to5*AJo  =  0-24,  X  =  3.55*  at  the  convective  wavenumber  and  the  primary 
energy  in  O(co)  results  from  turbulence  activity  in  the  outer  region  of  the  boundary  layer. 
However,  the  spectral  levels  at  wavenumbers  above  and  below  the  convective  ridge  are  directly 
influenced  by  the  small-scale  turbulence.  The  primary  energy  in  small-scale  turbulence  with  X  < 
lOOv/ux  exists  at  high  convective  wavenumbers,  kcv/ux  >  0.06,  and  corresponding  high 
frequencies,  cov/ux^  >  1-0.  Farabee  and  Casarella  concluded  that  energy  at  cov/ux^  >  0.3  results 
from  small-scale  turbulence  activity  in  the  buffer  region  of  the  boundary  layer.  Such  small-scale 
turbulence  has  little  effect  on  the  total  energy  in  the  autospectmm  O(co)  at  low  frequencies,  but 
significantly  influences  the  distribution  of  energy  in  the  wavenumber-frequency  spectrum  at  low 
frequencies.  This  conclusion  is  also  consistent  with  that  of  Farabee  (1986).  By  considering  the 
spectral  solution  to  the  Poisson  equation,  he  concluded  that  contributions  (of  various  levels)  to 
0(03)  at  low  frequencies  are  possible  from  turbulent  velocity  fluctuations  throughout  the  entire 
boundary  layer. 

Control  of  the  small-scale  turbulence  could  therefore  lead  to  improvements  in  flow-induced 
vibrations  and  noise  by  reducing  the  low-wavenumber  spectral  levels.  A  sharper  roll-off  for  small 
separations  was  shown  to  result  in  lower  spectral  levels  at  low  wavenumbers.  Control  aimed  at 
reducing  the  convection  velocities  of  the  small-scale  structures  could  therefore  lead  to  this  effect. 
However,  changes  in  the  spatial  decay  as  well  as  in  the  inherent  coupling  between  the  inner  and 
outer  regions  of  the  boundary  layer  will  clearly  have  an  impact.  Our  results  also  emphasize  the 
requirement  for  extremely  small  pressure  sensors  of  diameter  25v/ux  or  less  to  accurately  measure 
the  convection  velocities  at  very  small  spatial  separations. 

The  effect  of  the  convection  velocities  of  the  large-scale  structures  on  the  spectrum  was  also 
determined  by  varying  the  constant  value  of  Uc(X)  at  large  separations.  The  largest  value  of 
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Uc(X)  for  Farabee  and  Casarella's  low-frequency  data  was  0.87  at  X  =  106.  Use  of  the  function 
a(X-i-l)^,  a  constant  value  of  0.87  for  X  >  105  rather  than  0.85  for  X  >  80  had  a  negligible  effect 
on  the  spectrum.  Using  a  lower  value  of  0.80  for  X  >  29  had  a  very  small  effect  on  the  convective 
and  higher  wavenumbers.  However,  at  subconvective  and  lower  wavenumbers,  the  spectral  levels 
display  an  apparent  oscillation  as  shown  in  figure  6.  These  oscillations  thus  result  from  changes  in 
Uc(X)  for  X  >  29.  For  Farabee  and  Casarella's  boundary  layer,  X  =  29  corresponds  to  ^  = 
5,460v/ux,  or  4.75.  Structures  with  X  <  2,730v/ux  or  ^  <  2.35  have  a  negligible  influence  on 
Uc(X)  for  X  >  29.  The  convection  velocities  of  turbulence  with  X  >  2.65  in  the  outer  region 
therefore  also  influence  the  spectrum  at  subconvective  and  lower  wavenumbers.  For  comparison, 
we  also  considered  a  constant  value  of  0.75  for  Uc(X)  for  X  >  9.  This  leads  to  a  shift  of  the 
convective  ridge  to  higher  wavenumbers.  The  shape  of  the  spectrum  agrees  with  the  model  of 
equation  (6)  as  expected,  since  Uc(X)  is  constant  over  a  large  spatial  region. 

The  sensitivity  of  the  spectrum  to  the  decay  constant  y  was  determined  for  -0.165  <Y<  -0.110, 
with  Uc(X)  =  a(X-i-l)b.  At  the  convective  ridge,  the  levels  increase  as  lyl  decreases,  as  shown  in 
figure  7.  At  higher  wavenumbers,  the  levels  decrease  slightly  as  lyl  decreases.  At  subconvective 
and  low  wavenumbers,  the  levels  decrease  significantly  as  lyl  decreases,  with  oscillations  similar  to 
those  obtained  by  decreasing  the  constant  value  of  Uc  at  large  separations.  The  value  of  -0.145 
corresponds  to  the  data  of  Farabee  and  Casarella  at  n  =  0.24,  and  the  value  of  -0.125,  which  is 
accurate  at  higher  frequencies,  is  commonly  used  for  modeling  purposes.  Although  a  value  of 
-0. 1 10  does  not  support  the  data  at  this  frequency,  the  resulting  oscillations  in  the  spectral  levels  at 
subconvective  and  low  wavenumbers  are  interesting.  Whether  such  oscillations  are  physically 
realizable  or  simply  result  from  an  inaccurate  curve  fit  to  the  data  is  not  clear  at  the  present  time. 

Keith  and  Barclay  (1993)  investigated  the  effects  of  a  large  eddy  breakup  device  (LEBU)  on  the 
wall  pressure  cross  spectra.  The  LEBU  caused  no  measurable  changes  in  the  convection 
velocities,  but  led  to  a  significant  loss  of  coherence  and  greater  spatial  decay  at  low  frequencies. 
Reductions  in  the  autospectral  levels  also  occurred  at  low  frequencies.  Using  the  model  given  here 
by  equation  (4),  they  found  the  normalized  wavenumber-frequency  spectrum  had  reduced  levels  at 
the  convective  ridge  and  increased  levels  at  lower  and  higher  wavenumbers.  The  results  of  figure 
7  support  the  same  trends,  although  the  shapes  of  the  spectra  are  different  due  to  the  variation  of 
convection  velocity  with  separation.  Keith  (1989)  investigated  the  effect  of  riblets  on  the  wall 
pressure  cross  spectra.  The  riblets  produced  no  measurable  change  in  the  convection  velocities  or 
the  autospectra,  but  led  to  increased  coherence  levels,  implying  reduced  spatial  decay  at  low 
frequencies.  The  results  of  figure  7  indicate  that  riblets  may  therefore  lead  to  increased  levels  near 
the  convective  ridge  and  reduced  levels  at  subconvective  and  low  wavenumbers. 
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Recent  numerical  simulations  of  low  Reynolds  number  fully  developed  channel  flow  provide  a 
detailed  description  of  the  space-time  wall  pressure  field.  The  characteristics  of  the  streamwise 
wavenumber-frequency  spectra  computed  from  these  simulations  were  recently  discussed  by 
Chang  et  al.  (1994).  The  Direct  Numerical  Simulation  "DNS3"  (with  Rx  =  180)  and  Large  Eddy 
Simulation  "LES2"  (with  Rx  =  171)  results  of  Chang  et  al.  (1994)  for  frequency  Q  =  0.24  are 
compared  with  the  present  results  in  figure  8.  Here,  the  spectra  are  presented  in  normalized  form 
using  the  convective  wavenumber.  A  value  for  the  convective  wavenumber  was  determined  from 
the  peak  of  the  convective  ridge  in  each  case.  The  present  results,  where  the  variation  of 
convection  velocity  with  spatial  separation  is  taken  into  account,  show  improved  agreement  with 
both  the  DNS  and  LES  results  at  both  superconvective  and  subconvective  wavenumbers,  in 
comparison  to  the  constant  convection  velocity  (Corcos)  results.  The  Corcos  result  clearly 
overpredicts  the  levels  at  subconvective  and  low  wavenumbers.  A  value  of  -0. 145  for  the  decay 
constant  at  this  frequency  was  found  to  yield  the  best  agreement  with  the  numerical  results,  which 
is  consistent  with  the  cross-spectral  measurements  of  Farabee  and  Casarella  (1991)  at  this 
frequency. 


Figure  8.  Comparison  With  "DNS"  and  "LES"  Numerical  Results 
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ESTIMATED  TWO-WAVENUMBER-FREQUENCY  SPECTRUM  <D*(Ki,K2,0) 


Thus  far,  we  have  focused  on  the  wavenumber- frequency  spectrum  O  (Ki,Q),  which 
represents  an  integral  over  all  spanwise  wavenumbers  K2.  The  two- wavenumber-frequency 
spectrum  d>*(Ki,K2,n)  provides  the  distribution  of  energy  in  K2  as  well  as  in  Ki.  In  particular, 
an  estimation  of  0*(Ki,K2,0)  is  required  to  determine  the  effects  of  finite-sized  sensors  on 
measured  spectra.  A  two-dimensional  array  of  pressure  sensors  is  required  to  make  a  direct 
measurement  of  0*(Ki,K2,Q)  or  to  make  cross-spectral  measurements  from  which 
0*(Ki,K2,0)  may  be  estimated.  The  Corcos  model  for  the  full  cross  spectrum  may  be  expressed 
as 


G(X,Z,n)  =:  4>(£2)e(Vl£2XAJcll-i(n»Ucl))e(|JIOZAJc2l) ,  (8) 

where  (3  is  the  spanwise  decay  constant.  The  normalized  nondimensional  wavenumber-frequency 
spectrum  is  then  given  by 


0*(Ki,K2,i2)  =  j  J  e('yl^^X/Uci|-i(-Ki  -f-  Q/Ucl)X)e(P  I^AJr.?.h(K2Z))Hyrl7. , 

_oo  — oo 

(9) 

which  with  Ud  =  Uc2  =  ^^/Kc  (Ud  and  Uc2  are  assumed  to  be  independent  of  X  and  Z)  may  be 
evaluated  in  closed  form  to  obtain 

<D*(Ki,K2,i2)  =  4yP/((Kc2(Y2  +  (i  .  Ki/Kc)2))(p2  -h  (K2/Kc)2))  .  (10) 

The  spectra  obtained  from  equation  (10)  for  fixed  frequencies  Q  =  0.1,  1.0,  and  10.0  are  shown  in 
figure  9,  where  a  value  of  0.85  valid  for  large  separations  was  taken  for  Uc,  with  y  =  -0.125  and  P 
=  -0.70.  In  this  analysis,  we  have  not  attempted  to  take  into  account  the  dependence  of  Ud  and 
Uc2  or  y  and  P  on  frequency,  as  we  are  primarily  interested  in  determining  the  effects  of  the 
variation  of  convection  velocity  alone  on  the  spectra.  The  value  taken  for  p  was  based  upon  the 
measurements  of  Willmarth  and  Roos  (1965).  The  spectra  are  normalized  in  wavenumber  and 
spectral  level  using  the  convective  wavenumber,  which  was  defined  at  each  frequency  to  be  the 
streamwise  wavenumber  at  which  the  spectrum  0*(Ki,f2)  reached  a  global  maximum.  The 
spectra  are  symmetric  about  the  convective  peak  in  Ki  and  about  the  origin  in  K2.  With  increasing 
frequency,  the  spectral  peaks  decrease  in  magnitude  and  the  distribution  of  energy  becomes 
broader  in  wavenumber. 
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Equation  (9)  was  evaluated  with  Ud  =  a(X+l)t>  (where  a  and  b  are  as  previously  defined)  and 
Uc2  =  0.85  for  all  frequencies.  There  are  presently  insufficient  experimental  and  numerical  data  to 
determine  the  effectiveness  of  different  forms  for  Uc2  for  collapsing  cross-spectral  data.  The 
absence  of  a  mean  spanwise  convection  velocity  means  that  Uc2  represents  a  scaling  parameter 
rather  than  a  measurable  quantity,  as  in  the  case  of  Ud.  Here  we  consider  only  a  constant  value  of 
0.85  for  Uc2-  For  frequencies  Q  =  0.1  and  1.0,  the  convective  peak  displays  a  significant 
asymmetry  in  Ki,  over  a  broad  range  of  K2,  as  shown  in  figure  10.  This  effect  is  not  as  apparent 
at  the  highest  frequency,  Q  =  10.0.  The  convection  velocity  Ud  enters  only  through  the  similarity 
variable  QXAJcl-  At  higher  frequencies,  the  effect  of  the  variation  of  Ud  with  separation  on  the 
spectra  is  therefore  diminished  due  to  the  increasing  values  of  Q. 

CORRECTIONS  TO  AUTOSPECTRA  MEASURED  WITH  CIRCULAR  SENSORS 

Corcos  (1963)  derived  an  expression  for  the  ratio  of  the  measured  to  the  true  autospectrum 
<I)ni(w)/^t(“)  by  formulating  the  problem  in  the  spatial  domain.  The  spatial  formulation  leads  to 
the  attenuation  function 

<l)m(co)/<I)t((o)  =  J  /  e(Yl®^Ucl(^’T]>“)|-i"^ud(^,ri,co))e(plc0r|/uc2(^.'n>“)l) 

—OO  —CO 

T(^,ri)d^dri  >  (11) 

where  the  sensor  spatial  correlation  function  is  given  for  a  circular  sensor  of  radius  r  by 
T(tTl)  =  { 2/k  r2 }  {cos- 1  (e/2r)-(e/2r)(  1  -  (e/2r)2)  1/2} ,  e/2r  <  1  , 

T(^,ti)  =  0,  e/2r>l, 


and 


e  =  (^2+T^2)l/2.  (12) 

Here  we  have  considered  the  most  general  case,  allowing  uci  and  uc2  to  vary  with  Ti,  and  O).  In 
this  case,  equation  (11)  cannot  be  evaluated  in  closed  form  and  requires  numerical  integration.  In 
view  of  equation  (12),  the  limits  of  integration  for  equation  (11)  must  extend  only  over  the  region 
bounded  by  e/2r  <=  1.  For  this  reason,  the  spatial  formulation  is  computationally  less  intensive 
than  the  analogous  wavenumber  formulation. 
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Figure  9.  Computed  Wavenumber-Frequency  Spectrum  (l)a*(Ki,K2,fi)  for  ud  =  Uc2  -  (^-85 


(gp)  pAoq  uinj]03ds  pszg^uuojs^ 
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Figure  10.  Computed  Wavenumber-Frequency  Spectrum  <5b*(Kl,K2,Q)  for  Uci(^)  =  a(^-hl)b 

and  Uc2(^)  =  a(ri-i-I)b 


In  formulating  his  results,  Corcos  assumed  Uc  was  independent  of  spatial  separation,  but  could 
vary  with  frequency.  For  that  case,  with  ucl(^,Tl,co)  =  Uc2(^7'n>  =  uc(co)  and  application  of  a 

change  of  variables  in  equation  (11),  the  attenuation  Oni(w)/^t(®)  be  shown  to  depend  only 
upon  the  quantity  cor/uc(CD).  The  attenuations  may  therefore  be  computed  in  terms  of  the  variable 
cor/uc(co),  without  specifying  a  particular  form  for  uc(co).  Our  results  were  computed  with  a 
constant  value  of  -0. 125  for  y  and  -0.70  for  P  in  equation  (1 1).  We  therefore  have  not  attempted  to 
take  into  account  the  dependence  of  y  and  P  on  frequency.  The  slightly  higher  levels  in 
comparison  to  those  of  Corcos  (1963),  as  shown  in  figure  1 1,  are  thought  to  result  from  the  values 
for  y  and  P  used  by  Corcos  (not  stated)  and  possibly  the  computational  methods. 

We  are  primarily  interested  in  determining  if  the  variation  of  convection  velocity  with  separation 
has  a  significant  effect  on  the  attenuations  in  the  frequency  spectrum.  In  the  previous  section,  we 
determined  <l>*(Ki,K2,0)  at  three  discrete  frequencies  and  did  not  attempt  to  take  into  account  the 
variation  of  convection  velocity  with  frequency.  Here  we  estimate  the  attenuations  over  a 
continuous  range  of  frequencies.  We  therefore  formulate  equation  (11)  so  as  to  express  the 
attenuations  in  terms  of  0)r/uco(^o>^o>«>)  with  uco(^o>Tlo>a>)  an  arbitrarily  chosen  "reference" 
convection  velocity,  which  may  vary  with  frequency.  We  further  make  the  simplifying  assumption 
that  the  convection  velocities  are  separable  in  frequency  and  space,  such  that 


Figure  11.  Autospectrum  Attenuations  for  Constant  and  Spatially  Varying  Convection  Velocity 
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ucl(^.^.«»)  =  g(w)hl(^,ri) , 


(13) 


uc2(^,^.<a)  =  g(w)h2(^,Ti) . 

Although  there  is  insufficient  experimental  data  to  rigorously  justify  this  assumption,  the  results  of 
Abraham  (1994)  obtained  from  the  DNS  channel  flow  database  (see  Choi  and  Moin  (1990)) 
suggest  a  separable  form  for  ucl(^,co).  We  also  assume  the  same  frequency  dependence,  g(co),  ^ 

holds  for  both  functions. 

We  define  a  transformation  of  variables  given  by 

X  =  ^r  ,  z  =  Tj/r ,  (0*  =  tor/uco(^o.^o.®)  , 

such  that  equation  (11)  becomes 

oo  oo 

9[)ni(co)/<I>t(co)  =  j  j  e(Tl“*xh(xo,Zo)/h(x,z)|  -ico*xh(xo,Zo)/h(x,z)) 

_oo  — oo 

e((3lo)*zh(xo,Zo)/h(x,z)l)T=!=(x,z)dxdz .  (14) 

With  this  transformation,  T*(x,z)  does  not  depend  explicitly  upon  r,  and  for  any  h(x,z)  and 
h(xo,zo),  the  attenuation  depends  only  upon  the  quantity  co*.  Here,  the  attenuation  depends 
implicitly  on  the  frequency  dependence  of  the  convection  velocity,  g(co),  only  through  the  variable 
CO*.  We  arbitrarily  choose  h(xo,Zo)  =  h(0,0)  =  0.645Uo.  We  assume  that  ucl(^,ri,co)  =  ucl(^,co) 
and  uc2(^.'n»w)  =  Uc2('n»t^)  such  that  the  cross  spectrum  retains  a  separable  form,  with  hi(^)  = 
a(^  +  1)^  and  h2(Tl)  =  a(Tl  +  1)*^  and  with  a  and  b  as  previously  defined.  In  this  case,  there  Was  no 
measurable  change  in  the  attenuation  function  with  respect  to  the  former  results,  as  shown  in  figure 
11.  The  computed  attenuations  for  six  different  values  of  r/5*  collapse  onto  a  single  curve  when  f 

scaled  with  to*,  consistent  with  the  formulation  given  by  equation  (14).  The  change  in  slope  for 
cor/uco  ~  4  reflects  the  first  zero  in  the  wavenumber  response  of  the  circular  pressure  sensor.  t 

The  interpretation  of  the  effects  of  sensor  size  on  attenuations  in  the  autospectrum  is  aided  by 
formulating  the  problem  in  the  wavenumber  domain,  as  shown  by  Keith  and  Bennett  (1991). 

Here,  the  redistribution  of  energy  with  wavenumber,  as  shown  to  occur  in  figure  10,  is  insufficient 
to  significantly  affect  the  attenuations  in  the  autospectrum,  due  to  the  inherent  averaging  over  all 
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wavenumbers  resulting  from  the  sensor.  The  increases  and  decreases  in  levels  are  therefore 
effectively  canceled  by  the  wavenumber  response  of  the  sensor. 

CONCLUSIONS 

Comparison  of  convection  velocities  from  experimental  measurements  and  from  numerical 
simulations  covering  two  orders  of  magnitude  in  R0  supports  an  outer  variable  scaling  for 
collapsing  the  data.  The  early  work  of  Corcos  provided  a  model  for  the  wall  pressure  cross 
spectrum  in  terms  of  a  similarity  scaling  based  upon  a  convection  velocity  assumed  independent  of 
spatial  separation.  Recent  measurements  still  support  this  similarity  scaling,  with  the  dependence 
of  uc  on  spatial  separation  taken  into  account.  An  improvement  in  the  wavenumber-frequency 
spectrum  0*(Ki,Q)  estimated  from  the  cross  spectrum  is  obtained  by  including  this  effect  when 
computing  the  spatial  Fourier  transform,  as  evidenced  by  the  favorable  comparison  with  recent 
numerical  results.  The  redistribution  of  energy  in  wavenumber  resulting  from  the  spatial  variation 
of  the  convection  velocity  did  not  affect  the  attenuations  in  the  autospectra,  due  to  the  inherent 
averaging  over  all  wavenumbers  by  the  sensor. 

In  view  of  the  sensitivity  of  the  estimated  spectrum  to  the  convection  velocity  and  also  the  decay 
constants,  the  limited  cross-spectral  data  presently  available  preclude  formulating  a  model  for  the 
wavenumber-frequency  spectrum  valid  over  a  wide  range  of  frequencies  as  well  as  Reynolds 
numbers.  We  have  used  Farabee  and  Casarella's  measurements,  which  provide  the  most  extensive 
set  of  streamwise  cross- spectral  data.  For  applied  problems  of  flow-induced  vibration  and  noise, 
estimates  of  the  waU  pressure  wavenumber-frequency  spectrum  may  be  determined  numerically  by 
use  of  the  methods  presented  here,  taking  into  account  effects  related  to  the  particular  frequencies 
and  Reynolds  numbers  of  interest.  Although  we  have  used  convection  velocity  data  here,  the 
actual  phase  data  will  improve  the  accuracy. 
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